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Abstract 

New measurements by the NA49 experiment of the centrahty dependence of event-by-event fluc- 
tuations of the particle yield ratios (K+-|-K~)/(7r+-|-7r"), (p+p)/(vr++7r~), and (K+-|-K^)/(p+p) 
arc presented for Pb+Pb collisions at 158A GeV. The absolute values of the dynamical fluctu- 
ations of these ratios, quantified by the measure (Jdyn, increase by about a factor of two from 
central to semi-peripheral collisions. Multiplicity scaling scenarios are tested and found to apply 
for both the centrality and the previously pubhshed energy dependence of the (K+-|-K^)/(7r+-|-7r~) 
and (p+p)/(vr"''+7r~) ratio fluctuations. A description of the centrality and energy dependence of 
(K+-|-K~)/(p+p) ratio fluctuations by a common scaling prescription is not possible since there is 
a sign change in the energy dependence. 



1. INTRODUCTION 



The search for structures in the QCD phase diagram, hke the first order phase transition 
hne from hadronic to partonic degrees of freedom or the critical endpoint, has become one 

nn 

of the main activities in current and future high-energy heavy- ion experiments [1H3|. The 
experimental signatures for these structures are the subject of ongoing discussions. Lattice 
QCD calculations show that in the co-existence region of hadronic and partonic degrees 
of freedom and in the vicinity of the critical endpoint event-by-event fluctuations of, for 
example, the strangeness-to-entropy ratio increase significantly |4-l7|. Thus, a measurement 
of the energy dependence of a quantity sensitive to this ratio and an observation of a non- 
monotonic behavior may provide an indication of the location of the critical endpoint. 

The NA49 experiment at the CERN Super Proton Synchrotron (SPS) analyzed the en- 
ergy dependence of the ratio of inclusive K"*" and vr^ yields in central Pb + Pb collisions 
and observed a peak structure at beam energies around 30 - 40 A GeV js]. This motivated 
the analysis of event-by-event fluctuations of the (K++K^)/(7r"'"+7r~) (denoted K/tt) {q], 
(p+p)/(7r++7r^) (denoted p/vr) [o?] and (K++K~)/(p+p) (denoted K/p) ratios 10| as func- 
tion of the center-of-mass energy by means of the observable (Jdyn(see see Eqs. HfSl in section 
3.4), which measures the dynamical contribution to the fluctuations of the event-by-event 
particle ratios. The K/vr ratio fluctuations show a continuous increase towards lower colli- 
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sion energies, which is not reproduced by the UrQMD model pj|, but obtained qualitatively 
by HSD model calculations 12|. The p/vr ratio fluctuations as a function of the center- 
of-mass energy show negative values which indicate strong correlations. This observation 
is well reproduced by UrQMD model calculations and can be interpreted as the result of 
the production of nucleon resonances and their decays into pions and protons. The K/p 
ratio fluctuations exhibit a change of sign at ^ 30^4 GeV beam energy which is not well 
understood [10]. In view of the complex energy dependence of the fluctuations of the three 
particle ratios an additional study of their collision centrality dependence at the top SPS 
energy may help to clarify the interpretation. In particular, such an investigation may help 
to distinguish the contributions of the changing multiplicities and the genuine energy and 
collision volume dependence of the underlying correlations jl3 |. 

The STAR collaboration at the Relativistic Heavy Ion Collider (RHIC) also published 
results on particle-ratio fluctuations jl4| employing the observable z/dyn (Eq. [7] in section 



3.4). First results from a recent low energy scan in Au+Au collisions at ^/sNN = 7.7 GeV 
were presented at conferences [isl and show a different trend for the energy dependence of 
K/tt and K/p fluctuations when compared using the equivalence relation between z/dyn and 
o^dyn (see Eq. [7] in section 3.4). However, acceptance in both rapidity y and transverse 
momentum as well as the selection procedure of collision centrality differ. 

In this paper we address the dependence of event-by-event fluctuations of particle yield 
ratios on the centrality of Pb+Pb collisions in a fixed acceptance and at a beam energy 



of 158A GeV 



16|. In section 2 we describe the experimental equipment, in section 3 the 



analysis procedures. Section 4 presents the experimental results and compares to various 
proposed multiplicity scaling schemes. A summary section 5 closes the paper. 



2. THE NA49 EXPERIMENT 

NA49 is a fixed target experiment [17} at the CERN SPS. The trajectories of charged 
particles are reconstructed in four large volume Time Projection Chambers (TPCs). Two 
of them (VTPCs) are placed inside of two superconducting dipole magnets for momentum 
determination. Two main TPCs (MTPCs) are located downstream of the magnets on both 
sides of the beam. The performance of the MTPCs is tuned for high precision measurements 
of the specific energy loss dE/dx, which is the basis for particle identification employed in 



this analysis (see section 13.31 and 18|]). Except for the trigger and beam intensity the 



The Pb beam had 



experimental conditions in this analysis are the same as described in 
a typical intensity of 10^ ions/s. The minimum bias trigger was derived from a He-Cerenkov 
counter placed behind the target. Only interactions which reduced the beam charge and 
thus the signal seen by this detector by at least 10%, were accepted. The trigger cross section 
thus defined is 5.7 b out of a total inelastic cross section of 7.15 b. The resulting ensemble 
of 174 K events was divided into centrality classes according to the energy measured in the 
Veto Calorimeter (VCAL) located 26 m downstream from the target. 



3. DATA ANALYSIS 



3.1. Event, track selection and acceptance 

In order to reject backgound interactions, a valid fit of an event vertex was required and 
a cut around tlie known target position was applied. The contamination by background 
events remaining after cuts on vertex position and quality amounts to less than 5% for the 
most peripheral collisions and is negligible for near-central collisions. 

The useful acceptance for pions, kaons, and protons is constrained by the needs of particle 
identification. The separation power is highest for particles with large track lengths in the 
MTPCs which limits the analysis to the forward hemisphere in the center-of-mass frame. The 
coverage in the azimuthal angle is a function of center-of-mass rapidity y and transverse 
momentum px- The loose and tight sets of track cuts used in the present analysis are given in 
Table [B These are identical to those employed previously in NA49 analyses of fluctuations jo]. 
The acceptance is not only determined by the track selection cuts. In addition, only phase 
space bins are used for which the inclusive dE/dx distributions have more than 3000 entries. 



TABLE I: Loose and tight set of track cuts used in the analysis. 



Cut description 




Cut 




Loose 


Tight 


{dE/dx)/{dE/dx)MW 


< L8 


< L8 


Number of points in MTPC 


> 30 


> 30 


Number of points in VTPCl 




> 10 


Number of points in VTPC2 




> 10 


Fraction of potential points found in MTPC 


> 50% 


> 50% 


Number of entries required in phase space 


>3000 


>3000 


bin for fit of inclusive dE/dx distribution 






Cut in proton rapidity for 0.2 GeV/c 


y ^ ybeam ~ 


1 2/ < 2/beam " 1 


Track fitted to primary vertex 




yes 


impact parameter x-projection 




< 4 cm 


impact parameter y-projection 




< 0.5 cm 



The acceptance after all selection cuts is shown in Fig. [T] for central colli- 



sions. The range in pt varies slightly depending on the number of events in the 




FIG. 1: (Color online) Acceptance in center-of-mass rapidity y and transverse momentum px for 
the most central Pb+Pb collisions at 158A GeV. 



3.2. Collision centrality determination 

The determination of the centrality of the collisions is based on the energy of forward 
going projectile spectators as measured in the VCAL. The distribution of the VCAL energy 
-EvETO together with the division into 5% bins of the total inelastic cross section is shown 
in Fig. 12 

The energy resolution of the VCAL measurement is dominated by two effects: the in- 
trinsic energy resolution as given by the longitudinal sampling structure and by the non- 
uniformit y of light collection efficiency. The overall resolution of the calorimeter was shown 
to follow |17] : 



with E in units of GeV. 

The choice of 5% centrality bin size is motivated by the energy resolution of VCAL, the 
requirement to keep the reaction volume fluctuations at a minimum and the necessity to 
have sufficient statistics in each centrality bin. Volume fluctuations are relevant for ratios 
involving kaons, since their multiplicity does not strictly scale with the number of wounded 
nucleons iVw, or equivalently the reaction volume, in contrast to the multiplicity of pions 




FIG. 2: (Color online) Distribution of the total energy -Eveto of the projectile spectators deposited 
in the VCAL of NA49 in Pb + Pb collisions at 158^4 GeV beam energy. An event vertex cut (see 
text) was applied to remove background triggers. Vertical lines separate bins of 5% of the total 
inelastic cross section. 



and protons [19]. The influence of volume fluctuations on all particle-ratio fluctuations was 
studied by varying the centrality bin widths in the range from 3% - 20%. The results shown 
in Fig. |3]for the example of K/vr ratio fluctuations in the most central collisions led us to 
choose 5% wide centrality bins, the smallest bin size that leaves sufficient satisfies. For each 



bin the corresponding average number of wounded nucleons (A^w) was obtained 



Glauber model approach using a simulation with the VENUS event generator 20, 



rom the 
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FIG. 3: (Color online) Dependence of the cTdyn measure of K/vr ratio fluctuations for the most 
central collisions on the width of the centrality bin in Pb+Pb collisions at 158^ GeV. 



3.3. Particle identification by dE/dx 



The event-by-event measurement of particle ratios ideally implies track-by-track identifi- 
cation of the different particle types. The NA49 experiment provides energy loss measure- 
ments along the particle trajectories in the MTPCs with a resolution of approximately 4 % 
in the relativistic rise region for particle momenta p above 3 GeV/c. Since the separation of 
the dE/dx signals of pions, kaons, and protons at a given momentum is of the same order, 
track-by-track identification of particle types is not possible. Instead, we employ a statistical 
method, namely the Maximum Likelihood Method (MLM) to extract particle ratios from 
event- wise dE/dx distributions of negatively and of positively charged particles. 

In a first step energy loss distributions of all accepted tracks in the event ensemble were 
constructed in bins of p, px and 0. The binning details are shown in Table [ITl 

TABLE 11: Binning in phase space used for fitting tfie inclusive dE/dx distributions. Due to 
overlap of the distributions for different particle species around momenta of 3 GeV/c the first 3 
bins in total momentum were not used. 



Variable 


Range 


-^bins 


Bin size 


P 


1-120 GeV/c 


20 


logarithmic 


PT 


0-2 GeV/c 


10 


0.2 GeV/c 


</> 


0-27r 


8 


0.25-7r 


charge q 


L-1 


2 





The resulting inclusive specific energy loss distributions in each phase space bin were 
fitted with four Gaussian functions all having the same width for electrons, pions, kaons, 
and protons (and their antiparticles) . The values of the nine fit parameters (eight positions 
and one width) define Probability Density Functions (PDF) for each phase space bin and 
were stored in a look-up table for later use in the event-by-event fits. 

Using the PDFs one can calculate for each particle the four probabilities to be an 
electron (positron), a pion, a kaon or a proton (anti-proton). The sum of these probabilities 
weighted with coefficients 6^ become the factors in the likelihood function which depends 
on the coefficients 6^: 



N 

L{{e^}) = m2^aUq\p\pi,<P\idE/dxy) , (2) 

1=1 a 

where the index i runs over the particles of the event. The coefficients 6a are the relative 
yield fractions of each particle type in the event. The sum of the weights is constrained to 
unity: 



(3) 



By maximizing the likelihood function with respect to the relative yield fractions one 
obtains the best estimate of the different particle multiplicities in a given event. More 
details about the employed MLM can be found in 22|. 



3.4. Extraction of dynamical fluctuations 

The fluctuations of particle ratios in the event ensemble are defined as the ratio of the 
root of the variance ^yVa,T{A/ B) of the distribution of the event-wise particle yield ratio 
A/B to the mean {A/ B) of the same distribution: 



"° WB) ■ 

Defined in this way Udata will contain contributions from the finite number statistics, de- 
tector resolution, nonperfect particle identification and the genuine dynamical fluctuations. 
The ffist three contributions are considered as background. Since their contributions dom- 
inate the ratio fluctuation signal, their magnitudes have to be determined quantitatively. 
For an estimate of the statistical fluctuations and the detector resolution effects the event 
mixing method was applied. A new ensemble of artificial events was generated which contain 
particles from different real events, selected randomly such that in each artificial event no 
pair of particles originates from the same data event. In addition the multiplicity distribu- 
tion of the mixed events was constructed to be the same as the corresponding distribution 
of the real events. By this token dynamical fluctuations, which may be present in data, are 
absent in the sample of mixed events. The measure o"mix, evaluated according to Eq. H] for 
the mixed events, contains thus only the background fluctuations. Examples of distributions 



of the event- wise particle ratio for real and mixed events are shown in Fig. H] for central and 
semi-peripheral Pb+Pb collisions. 
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FIG. 4: (Color online) Event-by-event particle-ratio distributions for central (left) and semi- 
peripheral (right) in Pb-|-Pb collisions at 158^ GeV. The data points show real event and the 
histogram mixed event distributions. Loose track cuts (see Table |T]) were applied. 



We now define dynamical fluctuations (o"dyn) as the geometrical difference between the 



fluctuations measured in real and mixed events: 



Alternatively, the observable u^yn |23|, defined as: 



l^dyn — — t'stat 

_ Var(A) Var(5) Cov{A,B) (6) 

has been used to measure dynamical fluctuations of the particle ratio A/B. Here Ustat = 

1/{A) + 1/{B) is the contribution from finite number statistics. Assuming that detector 



effects cancel in adyn it was shown that adyn is related [13|, llJ] to the fluctuation measure 

^dyn- 



sign(adyn)a^yn ^ ^dyn = ' + ' ' ' -2 ^^^^^^ . (7) 

For a check of the sytematic uncertainties inherent in the mixed event background subrac- 
tion procedure we also determined v^yn from our data. Owing to our non-perfect particle 
identification we again use mixed events to account for the background in the evaluation of 
i^dyn from the event-by-event fitted particle multiplicities: 

t'dyn = t'data " t'mix • (8) 

The resulting values for u^y^ and a^y^ were found to satisfy the equahty of Eq. [7] within the 
systematic uncertainties estimated for adyn- 

As can be seen from Eqs. [5] and [7] the values of adyn and fdyn can be positive as well 
as negative. Assuming Poissonian single particle distributions, correlations lead to negative 
values of adyn, while positive values are indicative of anticorrelations between the particles. 



3.5. Systematic error estimation 

In order to study the systematic uncertainties introduced by the track selection the results 
from the tight and loose sets of cuts (see Table [T]) were compared. We take the absolute 
difference between the results of the analysis with the two extreme conditions as an estimate 
of the corresponding systematic error. 

Other sources of systematic uncertainty for the determination of the particle-ratio fluc- 
tuations are the dE/dx resolution and the method of event-by-event particle identification. 



This systematic effect was studied with the help of simulated events from the UrQMD 
model 111]. In a first step the particles from the generated events were filtered by an accep- 



tance table, representing the phase space bins of the real data which had sufficient statistics 
for successful fits of the dE/dx distribution. Then for each accepted track a dE/dx value 
was randomly generated from a parametrization of the inclusive dE/dx distribution for the 
true particle identity which depends on particle type and phase space bin. Finally, the ac- 
cepted tracks with simulated dE/dx values were processed by the same analysis routines 
as the tracks from real data. In addition to the determination of particle multiplicities by 
the MLM the true particle identities as generated by the Monte Carlo code were stored. 
Figure [5] presents a comparison of the values of the dynamical particle-ratio fluctuations as 
obtained by using Monte Carlo identity and results from the dE/dx fit. 
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FIG. 5: (Color online) Centrality dependence of the measure adyn for K/vr (a), p/vr (b) and K/p 
(c) ratio fluctuations evaluated for events simulated by the UrQMD model [11'] using either true 
identity or the event-by-event fit results based on the simulated particle dE/dx. 



The comparison of the results from both identification methods suggests the particle 
identification method based on the event-by-event MLM fit of the dE/dx distributions is 
valid in the (A^w) range above 190. The difference observed for the K/p ratio was included 
in the systematic error. 



4. EXPERIMENTAL RESULTS AND DISCUSSION 



4.1. Centrality dependence of dynamical particle-ratio fluctuations 

In this section we present our results on the centrahty dependence of dynamical fluctu- 
ations of K/tt, p/tt, and K/p ratios in Pb + Pb collisions at 158^4 GeV (numerical values 
are listed in Table IIIip . The dot symbols in Fig. show the dependence of (Xdyn (mean value 
of the results for tight and loose track cuts) of the three ratios on the average number of 
wounded nucleons (A^w)- The systematic errors are indicated by the shaded bands. Also 
shown by square symbols are the values of dynamical fluctuations in central Pb + Pb colli- 
sions at 158^4 GeV beam energy from previous NA49 analyses j^, [l^, which used a different 
event ensemble. The results from both analyses are in good agreement. 
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FIG. 6: (Color online) Centrality dependence of the measure (Tdyn of K/vr (a), p/vr (b), and K/p 
(c) ratio fluctuations. Dots show results from this analysis, squares show previously published 



measurements [lol. The curves depict predictions of the UrQMD model 
acceptance. The shaded bands show the systematic errors. 



nj] for the NA49 



One observes the same trend for all considered particle ratios, namely that the absolute 
value of the dynamical fluctuations increases with decreasing centrality (decreasing N^^). 
Interestingly, the UrQMD model [ll| reproduces this behavior for all three ratios as demon- 
strated by the lines in Fig. [61 The model was previously found to fail in describing the 
energy dependence of crdyn(K/7r) and (Jdyn(K/p) whereas it reproduced cxdyn 



TABLE III: Numerical results for crdyn(K/7r), 

'7dyn(p/'''")i and (7(jyn(K/p) wltli statistical and sytem- 
atic uncertainties for seven centrality intervals in Pb+Pb collisions at 158A GeV. Also listed are the 
corresponding average number of wounded nucleons (Nw) and the average numbers of identified 
particles (vr), (K), (p) in the acceptance used for the analysis. 





0-dyn(K/7r) 


f^dyn(p/7r) 


o-dyn(K/p) 


(7r+ ± TT^ 


(K+ ± K- 


> (P±P) 


384 


3.7 ± 0.8 ± 1.2 


-4.9 ±0.4 ±2.9 


-5.0 ±0.7 ±4.0 


349.4 


45.4 


51.1 


352 


4.0 ± 1.1 ± 1.2 


-5.7 ±0.5 ±2.9 


-7.0 ±0.7 ±4.0 


284.5 


35.3 


39.3 


319 


6.1 ±0.9 ±1.2 


-7.1 ±0.6 ±2.9 


-9.2 ±0.8 ±4.0 


234.0 


27.7 


31.1 


286 


3.0 ± 1.7 ± 1.2 


-6.4 ±0.8 ±2.9 


-9.0 ± 1.7 ±4.0 


191.7 


21.6 


23.1 


253 


7.0 ± 1.5 ± 1.2 


-8.6 ±0.8 ±2.9 


-13.5 ±0.9 ±4.0 


156.4 


16.5 


18.1 


220 


6.9 ± 1.9 ± 1.2 


-8.7 ± 1.1 ±2.9 


-17.6 ± 1.9 ±4.0 


124.1 


12.0 


13.7 


193 


11.5 ±1.5 ±1.2 


-7.6 ±2.7 ±2.9 


-16.9 ± 1.9 ±4.0 


97.7 


8.5 


11.3 



4.2. Scaling behaviour of dynamical fluctuations 

In this section we discuss various multiplicity scaling prescriptions which were pro- 



posed [13|, [ij, |2J] with the aim of separating effects of changing average particle multiplicities 
from the energy and collision volume dependence of genuine dynamical fluctuations. It is 
important to note that for comparisons of experimental data with scaling calculations the 
measured multiplicities inside the experimental acceptances should be used. The analysis 
will be applied simultaneously to the centrality dependence reported in this paper and the 



aiitvaf 



energy dependence previously pubhshed in 

In [l^ it was shown that adyn is expected to have a strong multiplicity dependence and 



might scale with ^yT/l^Aj~^^TJ{B) , where (A) and {B) are the average numbers of accepted 
particles of type A and B. As shown in Fig. [71^a,b) the measurements of the centrality and 
energy dependence of K/vr and p/vr fluctuations are consistent with the proposed scaling 
(sometimes called Poisson scaling). This result suggests that a large contribution to the 
observed variations appears to be caused by the changing multiplicities rather than by 
changes of the underlying correlations. 

In contrast, the energy and centrality dependence of K/p ratio fluctuations, plotted in 
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FIG. 7: (Color online) Dynamical fluctuations of the K/vr (a), p/vr (b), and K/p (c) ratio as a 
function of ^Jl/{K) + l/(7r), ^Jl/{-9) + l/(vr), and ^Jl/{K) + l/(p) respectively, (vr), (K), and 
(p) are the average number of kaons and protons in the acceptance. The solid lines show fits to 
Poisson multiplicity scaling a^yn oc ^Jl^A)~^\^Yf(B) (see text). Shaded bands indicate systematic 
uncertainties. 

Fig. Wis) a function of a/1/ (K) + l/(p), are not compatible with a common multiplicity 
scaling. The energy dependence shows a change of sign, indicating a change in the underlying 
correlation around SOA GeV beam energy. On the other hand, the centrahty dependence 
exhibits a smooth decrease which is close to the Poisson multiplicity scaling behaviour (solid 
line in Fig. El^c)). 

As already mentioned in the introduction section the STAR collaboration has presented 
results on the collision energy dependence of particle-ratio fluctuations measured at RHIC 
in Au+Au collisions in terms of the observable Vdyn- First results for ^/snn = 7.7 GeV [l5 | 
presented at conferences show a trend which differs from our results for K/vr and K/p when 
compared using the relation between z/dyn and a^yn (see Eq. [7] in section 3.4). Intensive 
disussion could not yet determine the cause of the difference. However, we note that the 
acceptance in both rapidity y and transverse momentum as well as selection of collision 
centrality are not the same. 

Another scaling behaviour of the dynamical fluctuations of the p/tt ratio was proposed in 
[l^ based on the hypothesis that these originate from the production and decay of nucleon 
resonances. Such decays introduce correlations between p and vr. Assuming that the variance 
terms in Eq. [7] can be neglected, the corresponding cidyn can be approximated by the following 



equation: 



^ Cov{A,B) miB})' 



with the parameter a expected to have the value 0.5. 




FIG. 8: (Color online) Dynamical fluctuations of the p/vr ratio as a function of the product of 
average numbers of protons and pions in the detector acceptance. 

The energy and centrahty dependences of the dynamical fluctuations of the p/vr ratio 
expressed as functions of the product (p)(7r) are plotted in Fig. [HI A fit of the data points 
with Eq. in] resulted in a parameters equal to a = 0.66 ± 0.12 for the energy dependence 
and a = 0.51 ± 0.03 for the centrality dependence. This experimental observation supports 
the hypothesis that the source of the p/vr ratio fluctuations is nucleon resonance production 
and decay. 

An alternative scaling hypothesis was also investigated for the K/vr ratio fluctuations. 
Since (K) ^ (tt), the dominating term in Eq. [7] for the dynamical fluctuations of the K/vr 
ratio may be the kaon variance term, provided the covariance term can be neglected. Fig. [9] 
shows the energy and centrality dependence of the K/vr ratio fluctuations versus the number 
of kaons (K) in the acceptance. The curves in Fig. [9] indicate that also the function 

/((K)) = a+^ (10) 

provides a good fit to both the centrality and energy dependence of K/vr ratio fluctuations 
with a = 2.4 ± 0.8 and b = 62.1 ± 16.6. 

In the most peripheral collision events fewer phase space bins are useable because of the 
lower multiplicities. We repeated the analysis by restricting the extraction of adyn to this 




FIG. 9: (Color online) Dynamical fluctuations of the K/vr ratio as a function of average number of 
kaons in the detector acceptance. Open square symbols show the results for the acceptance of the 
most peripheral set of events. The solid line shows the fit of the centrality dependence with the 
function of Eq. [TOl 

smaller acceptance for all centralities. The multiplicities (K) in the restricted acceptance, 
of course, decrease. Nevertheless, the open square symbols in Fig. [9] demonstrate that the 
results for a^yn still follow the scaling of Eq. [TOl 

Presently the NA49 collaboration is in the process of developing and applying a new anal- 



ysis procedure (identity method 25|) for the determination of event-by-event particle ratio 
fluctuations. It is designed to unfold the second moments of the multiplicity distributions 
of protons, kaons and pions. With this information, more direct tests of various models will 
become possible. 



5. SUMMARY 



We presented new measurements of the centrality dependence of p/vr, K/vr, and K/p 
particle ratio fluctuations in terms of o"dyn obtained by the NA49 experiment from Pb + 
Pb collisions at 158 A GeV. The measure cxdyn increases in absolute value with decreasing 
centrality for all these ratios. Comparisons to various multiplicity scaling schemes were 
made to both the centrality and the previously published energy dependences. Fluctuations 
of the p/vr and K/vr ratios are consistent with Poisson multiplicity scaling, thus suggesting 
that changing multiplicities rather than varying genuine correlations are the main source of 
these dependences. The p/vr ratio fluctuations also scale with l/((p)(vr))°-^, supporting the 



assumption that they originate from production and decay of nucleon resonances. The K/tt 
ratio fluctuations are also compatible with a 1/ (K) behavior, suggesting fluctuations of the 
kaon multiplicity as the main source of the measured energy and centrality dependences. 
In contrast, multiplicity scaling cannot describe the measurents of K/p fluctuations consis- 
tently. Although the centrality dependence of the absolute value of K/p ratio fluctuations 
exhibits a smooth increase for more peripheral collisions and is compatible with Poisson 
multiplicity scaling, a sign change is observed for the energy dependence. Therefore, the 
correlations causing K/p fluctuations appear to be changing in the SPS energy range. 
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